Silver transport in Gea;Sei_2;:Ag materials: ab initio simulation of 

a solid electrolyte 

De Nyago TaferQ and D.A. DraboJ 

Department of Physics and Astronomy, 
Ohio University, Athens OH 4 5701, USA 

M. Mitkova 
Center for Solid State Electronics Research, 
Arizona State University, Tempe, AZ 85287-6206, USA 
(Dated: February 2, 2008) 

Abstract 

In this paper, we present models of Ge-Se glasses heavily doped with Ag obtained from ah 
initio simulation and study the dynamics of the network with an emphasis on the motion of Ag+ 
ions. The models are analyzed with partial pair correlation functions, static structure factors and 
novel wavelet techniques. The electronic properties are characterized by the electronic density of 
states and analysis of specific electronic eigenstates. As Ag content increases, the optical band gap 
increases. Ag diffusion is observed directly from thermal simulation. The most diffusive Ag"*" ions 
move preferentially through low density regions of the network and the existence of well-defined 
trapping centers is confirmed. Preliminary information about temperature dependence of trapping 
and release is provided. 

PACS numbers: 61.43.Fs, 66.30.Dn, 66.30.Hs, 71.23.An, 71.23.Cq, 71.15.Mb 
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I. INTRODUCTION 



The chalcogenide glasses are preferred semiconducting materials for applications. They 
have well-defined niches in fiber optics^, optical recording^, phase change memory^, and 
other technologies. Ge-Se glasses have been particularly studied because of their ready glass 
formation, easy synthesis requirements, and good chemical stability. The basic structural 
units are Se chains and Ge-Se tetrahedra which may combine in a variety of ways. Defects 
(including homopolar bonds) exist in these systems as seen in both experiment and theory. 

The GeSe system was the first in which formation of an intermediate phase was demon- 
strated experimentally by Boolchand et and further developed theoretically by Thorpe, 
Jacobs, Chubnysky and Phillips^. In binary Ge^jSci-j; glasses, the self-organized phase exists 
in the 0.20 < a: < 0.254 range, with glasses at a: < 0.20 regarded as floppy while those with 
X > 0.26 stressed rigid. Dynamic calorimetry measurements on the intermediate phase have 
led to the conclusion that such materials do not age&, a feature that may be of importance 
in application of these materials. 

Silver added to chalcogenide glass hosts has attracted widespread interest in soft con- 
densed matter science^. The interest emerges in part from the extensive bulk glass forming 
tendency in the Ge-Se-Ag ternary, the spectacular enhancement (eight orders of magni- 
tude) in electrical conductivity of glasses with Ag relative to the glassy chalcogenide hosts, 
and from light-induced effects such as photo-doping, photo-diffusion and photo-deposition. 
Although the mobile ions in amorphous materials have been studiect^ their detailed dynam- 
ics in amorphous hosts still constitutes one of the unsolved problems of solid state ionics. 
The structure of Ge-Se-Ag glass has been investigated using several experimental methods, 
including X-ray diffractio n^^i^^ neutron diffraction with isotopic substitution'^^, EXAFS^, 
differential anomalous X-ray scattering (DAS)^^*i^»i^ and Modulated Differential Scanning 
Calorimetry (MDSC) and Raman spectroscopy^!^. Despite this impressive database, the 
structure of the ternary Ge-Se-Ag glasses has not yet been completely determined. There 
continues to be a debate on basic aspects of the glass structure (i.e. homogeneity and Ag 
coordination) especially for Se rich glasses with more than 67% Se. 

Experimental evidence for macroscopic phase separation in these materials has come from 
MDSC results, which indicate bimodal glass transition temperatures^ . In these experiments, 
one Tg is independent of glass composition, and identified with a Ag2Se glass phased, while 
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the second that varies with glass composition is related to the Ge-Se backbone. For the 
time scales (and possibly also lengths scales) of our simulations, these effects do not emerge, 
but are an interesting challenge for the future. 

In a recent Letter^^, we briefly reported the motion of Ag ions in glassy chalcogenide hosts, 
and demonstrated the existence of ion trapping centers, which are important for relaxation 
processes in disordered systems^^. Here, we provide detailed information about structural 
and electronic properties, and also give new information about temperature dependence of 
the trapping, and also the geometry of the traps. In this paper, we have focused upon 
Ag-doped glasses containing Ge 25 at.% and Se 75 at.% (we later call this GeSes). This 
composition is near the intermediate phase (slightly into the stressed-rigid phase). To our 
knowledge, this is the first ab initio simulation of these materials. 

The rest of this paper is organized as follows. In Section II, we describe the simulation 
procedure to fabricate the atomistic models, discuss the ab initio total energy functional 
and force code used and other approximations. In Section III, we describe the structural 
properties using conventional measures such as static structure factors, and also apply a 
novel wavelet method to explore intermediate range order. Electronic properties are briefly 
discussed in Section IV, and Section V is concerned with the dynamics of the Ag"*" ions in 
the amorphous matrix. 

II. MODEL GENERATION 

A. Energy Functional and Interatomic Forces 

For the simulations reported in this paper, we use FIREBALL2000 developed by Lewis 
and coworkera^ii. This code is an approximate ab initio density functional approach to the 
electronic structure, total energies and forces based upon pseudopotentials and a real-space 
local basis of slighty-excited pseudoatomic orbitals to represent the Kohn-Sham functions. 
The method uses separable pseudopotentials, and allows the use of double-zeta numerical 
basis sets and polarization orbitals. The calculation is undertaken entirely in real space, 
which provides substantial computational efficiency. The exchange-correlation energy was 
treated within the LDA, using the results of Ceperley and Alder^- , as interpolated by Perdew 
and Zungei»2^ (more intricate gradient corrected functionals are available if needed). The 
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pseudopotential and pseudoatomic wave functions were generated in the Troullier-Martins 
forn>2^ employing the scheme of Fuchs and Scheffer-^i. Hamihonian and overlap matrix ele- 
ments are precalculated on a numerical grid and the specific values needed for a particular 
instantaneous conformation are extracted from the tabulated values via interpolation. Nat- 
urally, the integral tables need to be generated only once, for a given set of atomic species, 
rather than performing quadratures "on the fly" during a MD run. 

B. Model Formation 

The models described here were generated using the melt quenching method. We be- 
gan by randomly placing atoms in a cubic supercell according to the desired composition 
[for (GeSe3)o.9oAgo.io 54 germanium atoms, 162 selenium atoms and 24 silver atoms; for 
(GeSe3)o.85Ago.i5 51 germanium atoms, 153 selenium atoms and 36 silver atoms] with the 
minimum acceptable distance between atoms 2 A. The size of the cubic cells was chosen 
to make the density of these glasses close to experimental data. The box size of the 240 
atom supercell of (GeSe3)o.9oAgo.io and (GeSe3)o.85Ago.i5 are respectively 18.601 Aand 18.656 
Awith corresponding densityl? 4.98 g/cm^ and 5.03 g/cm^. The structures were annealed 
and we obtained well thermalized melts at 4800K. We took three steps to cool the cells. 
First, the cells were equilibrated to 1100 K for 3 ps; then they were slowly cooled to 300 K 
over approximately 5 ps. The MD time step was 2.5fs. Simple velocity rescaling was used for 
the dissipative dynamics. In the final step, the cells were steepest descent quenched to OK 
and maximum forces smaller in magnitude than 0.02 eV/A. All calculations were performed 
at constant volume using the F point to sample the Brillouin zone to compute energies and 
forces. 

The use of MD to quench a model liquid is natural and in principle completely gen- 
eral, since it superficially mimics the process of glass formation. It has been utilized for 
structural calculations for GeSe^^, As2Se3^, Ge-Se-Ag^'' and other systems (though there 
are sometimes hints that there is too much "liquid-like" character frozen into the resulting 
models). With accurate force calculations, simulated melt quenching has produced disap- 
pointing results for liquids such as GeSe^^ and ternary glasses such as As-Ge-Se^. We 
believe that the success with (GeSe3)i_2:Ag2^^ (x=0.10, 0.15) glasses is connected to the fact 
we are in a weakly overconstrained glass forming part of the vibrational phase diagrant^. 
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TABLE I: Basic short range order parameters in models (5^(GeSe3)o.9oAgo.io (2;=0.10) and g- 
(GeSe3)o.85Ago.i5 (2;=0.15) and compared to available data, ri is the average bond length, n the 
average coordination number, and r2 the second nearest neighbor distance. 







This Work 






Experiment (Ref»i^) 


Glass ri 


(±0.01 A) 


r2 (±0.01 A) 


n (±0.01) n 


(±0.02 A) 


ra (±0.05 A) n 


x=0.10 


2.47 


3.78 


2.525 


2.46 


3.80 2.471±0.002'^ 2.543±0.100 


x=0.15 


2.48 


3.82 


2.775 


2.48 


3.82 2.855±0.002" 2.852±0.100 



''Analytical calculation 
^Monte Carlo method 



The "cook and quench" method is less effective in the highly overconstrained regime. It is 
hkely that hybrid schemes mixing experimental information and ab intio simulation (such 
as the "ECMR" method^i) would be most effective in this composition regime. 

III. STRUCTURAL PROPERTIES 
A. Short range order and defects 

We define short range order (SRO) as a length scale not longer than the second nearest 
neighbor distance. Table HI gives an overview of the SRO in g'-(GeSe3)o.9oAgo.io and g- 
(GeSe3)o.85Ago.i5. The average bond length and average coordination number closely agree 
with the available datai^. With increasing Ag concentration, the average coordination num- 
ber increases to 2.8 in g'-(GeSe3)o.85Ago.i5 from 2.51 in g-GeSe^. Also listed is the second 
nearest neighbor distance. Table HTl lists the average bonding distances of different possible 
bonds present in the model. Different values of the Ag-Ag bond distance have been proposed. 
Using the differential anomalous scattering (DAS), Westwood et al^'^ obtained a value of 3.35 
A for the Ag-Ag distance, a bit longer than our observation or other data^^. Both models 
contain structural defects. In addition to the normally coordinated Ge4 and Sea: Ge3, Ses 
and Sei are present in both models. Table HTTl summarizes the statistical distribution of the 
main structural components. A look at the table shows that Ge-Se, Se-Se, Ag-Se and Ag-Ge 
correlations depend upon x with additional Ag modifying the Ge-Se and Se-Se bonding. By 
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integrating the partial pair correlation function ga/3(r) (Fig. HJ, we estimate the average 
coordination number of Ag in the models. We found an average coordination of 2.0 and 
2.9 in 5'-(GeSe3)o.9oAgo.io and 5'-(GeSe3)o.85Ago.i5 respectively. We note that this is hard to 
estimate without ambiguity because of the lack of a well-defined deep minimum in g{r) after 
the first peak especially for the x = 0.15 model. 

When discussing the coordination of Ag one can follow the arguments of Kastner— who 
was the first to describe the bonding of metal atoms in chalcogenide glasses. However, we 
have to consider the fact that in addition to the covalent bond that is expected to form 
with the chalcogens, Ag offers three empty s —p orbitals and is surrounded by the lone-pair 
electrons of its chalcogen neighbors. The latter offer the opportunity for the formation of up 
to three coordinate bonds. A coordinate bond is similar to a covalent bond and has similar 
strength but the bonding electrons are supplied by one bonding partner (the chalcogen 
atom)^^. As a result, the lowest energy-bonding configuration for Ag in chalcogenide glasses 
is an overall neutral complex with Ag positively charged with the negative charge located 
on neighboring chalcogen or chalcogena^. The coordination of Ag thereby can vary, but 
the expected value would be in average close to three - with one covalent bond and up to 
2 coordinated bonds. The opportunity for four-fold coordination also exists as there is one 
more free s — p orbital at the Ag atom but evidently this does not satisfy the requirements 
of the lowest energy configuration and the electronegativity of the entire complex in the 
presence of other cations in the system so that the probability for this type of bonding is 
lower. Our simulations appear to be consistent with these chemical considerations, as all 
of the Ag is two-fold in the 10% model, with the addition of three-fold Ag in the 15% Ag 
model. 

B. Intermediate range order 

Fig. [T] shows the calculated static structure factors for (GeSe3)o.9oAgo.io and 
(GeSe3)o.85Ago.i5 and the comparison with the experimental data from ref.— . As we did 
not include a priori information in the model formation process, the fact that the peak 
positions and spectral weight of S(Q) agree well with experimental data is encouraging. The 
third and fourth peaks are a result of the short range order in the models. The position of 
the second peak does not depend strongly on the Ag concentration (though its width does). 
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This peak is located at Q=2.09 A ^ in both models. By contrast the third peak inten- 
sity decreases with increasing Ag. This peak is located at Q=3.49 A^^ in (GeSe3)o.9oAgo.io 
and at Q^3.41 A'^ in (GeSe3)o.85A go.is; whereas in GeSca {x=0, concentration of Ag in 
(GeSe3)i_a;Ag2.) this peak appears near^^*^ 3.53 A^^. In Fig. Elthe partial structure factors 
show that SGeGe(Q) and SceSelQ) tend to cancel each other. 

Both models reveal a feature in S(Q) near 1.07 A"^. This is a harbinger of a First Sharp 
Diffraction Peak (FSDP), that becomes explicit in more Ge-rich materials. Piarristeguy et 
aU^ show that this peak varies as a function of Ag content. As Ag concentration increases, 
the FSDP decreases due to a change of the IRO. Moreover, Ag disturbs the GeSe4/2 network, 
and leads to the fragmentation of GeSe4/2 tetrahedrons. From the partial structure factors 
it is apparent that the "proto-FSDP" has contributions from all of the partials. To elucidate 
the intermediate range order more quantitatively, we have used the wavelet-based methods 
of Harrop and coworkers^. This is a promising scheme for interpreting structure factor data 
based upon continuous wavelet transforma^l (CWT). In Fig. El we illustrate the results of 
the wavelet analysis extracted from the experimental data^^. The most obvious feature is 
the extended range real-space correlations associated with the diffraction peak near 3.5A^^ 
for the 10% Ag model. This is connected to the narrowness of the peak for the 10% glass 
relative to the 15% material: in the latter case the correlations disappear by about 15A~^, 
whereas the correlations extend to at least 25Afor the 10% glass. A similar state of affairs 
also accrues for the peak near 5.5A~^ and for a similar reason. This work emphasizes that 
simple associations of the reciprocal of the peak position of the FSDP to real-space length 
scales is misleading, particularly if the peak is narrow as discussed by Uchino and coworkers 
in studies of silica glasses,. 

In Fig. |3]we plot the partial pair correlation function of both models. The Ge-Se and Se- 
Se pairs provide the dominant contribution to the first shell of the pair correlation function 
g(r) whereas Ag-Se contribute to the second peak, Ag-Ag to the third and Se-Se (second 
nearest neighbors) to the fourth peak. 

IV. ELECTRONIC PROPERTIES 

Having studied structural properties, we now briefly analyze the electronic properties of 
our models. The electronic density of states (EDOS) of both models are calculated and 
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TABLE II: Nearest-neighbor distances in (GeSe3)o.9oAgo.io (a;=0.10) and (GeSe3)o.85Ago.i5 
(x=0.15) glasses. 

Correlation Distance (±0.03 A) Ref.J^ (±0.05 A) 

a;=0.10 Ge-Ge 2.38 

Ge-Se 2.37 2.37 

Ge-Ag 2.35 

Se-Se 2.36 2.37 

Ag-Se 2.63 2.67 

Ag-Ag 3.10 3.05 

x=0.15 Ge-Ge 2.38 

Ge-Se 2.37 2.37 

Ge-Ag 2.35 

Se-Se 2.36 2.37 

Ag-Se 2.63 2.67 

Ag-Ag 3.21 3.05 



TABLE III: The statistical distribution of the main structural components in 5-(GeSe3)o.9oAgo.io 
(x=0.10) and g-(GeSe3)o.85Ago.i5 (x=0.15) models. Percentage of a given component in the total 
configuration is given in parentheses. Here the subscript number indicates the coordination number. 
For example, Ag2 means the two-bonded Ag sites. 





Ge4 


Gea 


Se2 


Se3 


Sei 


Ag2 Ag3 


x=0.10 
x=0.15 


35 (65%) 
34 (66.7%) 


17 (31.5%) 
13 (25.5%) 


86 (53.1%) 
85 (55.6%) 


56 (33.9%) 

57 (37.2%) 


17 (11.8%) 
11 (7.2%) 


24 (100%) 

31 (86.1%) 4 (11.1%) 




Ge-Se 


Se-Se 


Ag-Se 


Ag-Ge 






x=0.10 
x=0.15 


64% 
57.1% 


21% 
19.2% 


13% 
18.2% 


1.5% 
5.2% 
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analyzed by the inverse participation ratio (IPR), which we denote by T. The EDOS are 
obtained by summing suitably broadened Gaussians centered at each eigenvalue. The IPR 

n=l 

determines the spatial localization of electronic eigenvalues. Here is the number of atoms 
in the model and q{n, E) is the MuUiken charge localized on atomic site n in a certain 
eigenstate E. Hence, X is a measure of the inverse number of sites involved in the state 
with energy E. For a uniformly extended state, the MuUiken charge contribution per site is 
uniform and X(E)=1/N. For an ideally localized state, only one atomic site contributes all 
the charge and X(E) = 1. Therefore a larger value of X means that the eigenstate is more 
localized in real space. 

In Fig. we report the EDOS and the species-projected density of states of 
(GeSe3)o.9oAgo.io and (GeSe3)o.85Ago.i5 glasses (our electronic eigenvalues have been shifted 
in order to place the valence band edge eigenvalue at zero). It should be noted that the 
spectra of both models are similar and closely related to the EDOS of Ge^Sei-^: (x > 0.15). 
With the addition of Ag into (/-GeSea, an intense peak, due to the Ag Ad electrons appears 
at about -3.47 eV as shown in Fig. |3| The valence band exhibits three features. The two 
lowest bands between -14.8 eV and -7.0 eV originate from the atomic 4s-like states of Ge and 
Se partially hybridized to form bonding states to Ag atoms. The next band lying between 
-7.0 and 0.0 eV contains p like bonding states of Ge and Se and d like bonding states of Ag. 
The peak in the topmost valence region is due to the lone-pair 4p electrons of Se atoms. The 
r point optical gaps of 5'-(GeSe3)o.9oAgo.io and g'-(GeSe3)o.85Ago.i5 are respectively of the 
order of 1.20 and 1.26 eV^. As the Ag content increases, the optical band gap increases. To 
our knowledge, experimental information about the EDOS of both systems is unavailable, 
so the curve in Fig. |S1 is actually a prediction. This is an interesting contrast to the work of 
Simdyankin and coworkera^^, who show that for low concentrations of Cu, the gap decreases 
with addition of Cu in AsS and AsSe glasses. Care is needed in comparing these results 
since the hosts and transition metals are different, and our models have far higher metal 
content. 

In order to connect localized eigenstates to particular topological/chemical regularities we 
plot in Fig. inithe IPR in the band gap region. We found that the localization in the valence 
band is mainly due to Se atoms and in the conduction band to Ge atoms in both models. A 
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close look at the localized states at the band edges shows that the localized states at the top 
of the valence band of g'-(GeSe3)o.9oAgo.io are mostly associated with two- and three-fold 
coordinated Se atoms with homopolar bonds, whereas the localization at the conduction 
band edge arises from overcoordinated Se associated with homopolar bonds and four-fold 
coordinated Ge connected to Se atoms involved with Se-Se homopolar bonds. By contrast, 
the top of the valence band of 5'-(GeSe3)o.85Ago.i5 is quite extended; the conduction band 
edge shows very few localized states due to the overcoordinated Ge and Se. This explains 
the results of Kawasaki^ showing dominance of the ionic conductivity related to Ag+ ions 
at these particular compositions. 

V. THE DYNAMICS OF SILVER IONS 

An outstanding feature of the materials is the high mobility of the silver in the com- 
plex host. The chemical explanation for the high diffusivity of Ag is its high quadrupolar 
deformability^. To "mine" information from our MD simulations, we begin by computing 
the mean-square displacement (MSD) functions for all of the atomic constituents: 

(^'W)« = TrE(|ri(t)-ri(0)r), (1) 

where () is an average over MD simulation time and the sums are over particular atomic 
species, a. The MSD were calculated for Ge, Se and Ag ions for both models and are shown 
in Fig. [7| The MSD of Ag ions increases rapidly with time, whereas that of the Ge and Se 
show a very slight slope. We chose a temperature of 1000 K to illustrate the diffusion, and 
later discuss behavior at selected lower temperatures. 

To explain the mechanism of diffusion of silver, we examine the trajectories of these 
particles. We obtained 2.5x10^* steps of time development, for a total time of 62.5 ps and a 
fixed temperature of 1000 K. Fig. |H1 illustrates 2D projections of trajectories of the most and 
least mobile Ag atoms in the x = 0.15 model. We notice that for short times, the MSD of 
the most mobile atoms increases due to the diffusive motion of Ag. At intermediate times, 
the atoms may be trapped in a cage formed by their neighbors, and at the longest times we 
can explore, they can escape such traps and diffuse again. Thus, our trajectories can largely 
be separated into vibration around stable trapping sites and hops between such sites. In 
both glasses, a fraction of Ag atoms move large distances (see Fig. In (GeSe3)o.9oAgo.io, 
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about 91.67 % of silver atoms move an average distance greater than 2.5 Afor a time scale of 
39 ps. Among them 25 % have an average displacement greater than 5 A. By contrast only 
4.2 % of Ag atoms move less than 2 A. On the other hand, about 89 % of Ag atoms move 
on an average distance greater than 2.5 Ain (GeSe3)o.85Ago.i5 for the same period of time. 
36.1 % of those atoms have an average displacement greater than 5 A. The most mobile 
Ag atoms move on an average distance of 8.2 A. These numbers illustrate the high ionic 
mobility of Ag ions in these complex glasses and are suggestive for a significant contribution 
of correlated hops of Ag"*" in the diffusion process. 

Based on these trajectories, thermal transport coefficients such as diffusion coefficients can 
be evaluated. Typically one uses either the Green-Kubo formula^^ where the time-dependent 
velocity autocorrelation (VAC) is integrated or the Einstein relation^ is employed, and the 
MSD is differentiated with respect to time. Since transport coefficients are equilibrium 
properties, the system must be properly thermalized before the transport properties can be 
estimated. Our simulations show that at time t>4ps the systems are well equilibrated . The 
Einstein relation for self-diffusion reads: 

(|r(t)-r(0)|^)=6Dt + C (2) 

where C is a constant, D is the self- diffusion coefficient and (|r(t) — r(0)|^) is the mean- 
square distance from initial position at time t, averaged over atoms of a given species. Direct 
simulation of the atomic trajectory and simple fitting yields C and D, and, in particular, 
estimates for the self-diffusion coefficient of Ag, D^g- The estimated values of D^g as a 
function of temperature are listed in Table IIVI These results are qualitatively reasonable 
when compared to the recent (room temperature) experiments of Urena et al'^'^ with the 
appropriate exponential activation factor included. The probability of correlated motion of 
Ag+ will of course increase with increasing Ag concentration in agreement with decrease of 
the activation energy for conductivity^. 

The dynamics is sensitive to the temperature. We performed additional simulations at 
temperatures ranging from 640K to lOOOK. In Fig. ^Jwe illustrate the hopping, which is 
qualitatively like the high temperature hopping, and the traps are very well defined. We 
note that even the most mobile Ag ions spend a substantial time in the traps, and appear 
to hop very efficiently (quickly) between traps. 

To further study the trapping centers'^^ we also obtained estimates of trap sizes and trap 
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TABLE IV: Estimates for the trap size r^^, trap lifetime i^^, and the self-diffusion coefficient Dj^g 
as a function of temperature T. 





T (K) 


rtr (A) 

LI \ / 


Ur (ps) 

LI \Jr / 


{cm? Is) 


X = 0.10 


640 


0.81 


5.89 






800 


1.41 


3.31 


1.29x10"^^ 




1000 


1.78 


2.73 


1.59x10-5 


X = 0.15 


640 


0.94 


5.58 


1.52x10-*^ 




700 


1.34 


4.67 


1.82x10-6 




800 


1.40 


3.80 


6.50x10-6 




1000 


1.84 


2.66 


2.06x10-5 



lifetimes as a function of temperature. In our calculations, we consider only particles that 
experience more than one traps. To calculate the size of the traps we enclose the particle 
trajectories of each trap in a sphere of radius rtr centered on the average position of the 
particle in the trap. Then, we determine the displacement of the particle trajectories in the 
trap with respect to the average trajectory, and average over all the displacements. The 
trap size rtr is then obtained by averaging over different traps. Knowing the trap sizes, 
their lifetimes can be easily determined. In Table IIV[ we give an estimate of the trap sizes 
and trap lifetimes as a function of temperature. As T increases, the trap sizes increase and 
their lifetimes decrease. The averaged trap radii range from about 0.8 to more than 1.8 A, 
and the averaged lifetimes extend from about 2.5 to more than 5 ps. The individual trap 
radii can be as large as 2.4 Aand the lifetimes as long as 7 ps. We have also used the MSD 
method^ to estimate trap sizes, and obtain results within a factor of ~ 2 from our simple 
geometrical approach. 

Further insight into the mechanism of diffusion in (Gea;Sei_2^)i_yAgj, can be found by 
studying the behavior of the molar volume of particular regions containing silver atoms. 
Hence we calculate the local density of the most and least mobile silver atoms as a function 
of time, then compare them to the density of the glass. To do so, we draw a sphere of radius 
R=4 A. The center of the sphere is the position of the Ag atoms we are tracking at a time t 
(the center of the sphere varies as a function of time). Then we calculate the mean density 
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of atoms inside the sphere. Fig. illustrates the local density of a few Ag atoms as a 
function of time. As seen on the figure, the most mobile Ag atoms are consistently located 
in regions with a lower local density (lower local volume fraction) and higher disorder. On 
the other hand, as we showed by direct calculation, there is little if any correlation between 
the trajectory-averaged mean density and the tendency to diffuse for the Ag atoms. On 
the other hand, perhaps unsurprisingly, we found also explored correlations between the 
average displacement of mobile Ag ions and the standard deviation of their local density 
{<Ji=\J (pf) — {piy, where pi is the local density of the Ag ion i, and () means trajectory 
average.). We found that as the average displacement increases, the standard deviation 
becomes larger (see Fig. IT^ . The correlations are perhaps linear with much noise. This 
implies that the diffusive Ag ions are exploring a wide variety of densities and the weakly- 
diffusing Ag sample a restricted density range. 

VI. CONCLUSION 

We have presented ab initio models of GeSe glasses heavily doped with Ag and studied the 
dynamics of the network with an emphasis on the motion of Ag ions. The models reproduce 
structural data, including reasonably subtle features in the diffraction data including the 
first peak (or shoulder) in S(Q). Wavelet methods help significantly in revealing intermediate 
range real-space correlations in the glasses. The atomistic motion of Ag+ ions is detailed for 
short times with a reliable first principles interaction. We have shown by direct calculation 
that trapping centers exist, and have shown that local basis ab initio MD can provide direct 
insight into the processes of transition metal dynamics in amorphous chalcogenide materials. 
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FIG. 1: Calculated total structure factor S(Q) of (GeSesjo.goAgo.io and (GeSe3)o.85Ago.i5 glasses 
compared to experimentii^. 



S. I. Simdyankin, M. Elstner, T. A. Niehaus, Th. Frauenheim and S. R. Elliott, 
cond-mat/050 4207' (2005) . 

The optical gap is underestimated by nearly a factor of 1 /2. This is an expected artifact of the 
Kohn-Sham eigenvalues. 

M. Kawasaki, J. Kawamura, Y. Nakamura, and M. Aniya, Sol. State Ion. 123, 259 (1999). 
David Chandler, Introduction to Modern Statistical Mechanics (Oxford University Press, New 
York, 1987), pp.246-252. 

M.A Urena, A. A. Piarristeguy, M. Fontana, B. Arcondo, Sol. State Ion. 176, 505 (2005). Using 
the diffusion coefficient of Ag at room temperature and the parameters obtained by these 
authors, we extrapolate the diffusion coeffient of Ag at 1000 K, 800 K, 700 K, and 640 K. We 
obtained diffusion coefficient of the order 1/10 times the theoretical values for T=1000 and 640 
K. In view of the short simulation times and small models, this is reasonably consistent with 
experiment. 

M. Gutenev, A. Tabolin, and T. Rykova, Phys. Chem. Glasses 17, 36 (1991). 

J. C. Phillips, Rep. Prog. Phys. 59, 1133 (1996). 

E. R. Weeks and D. A. Weitz, Phys. Rev. Lett. 89, 095704 (2002). 



from Flying Frog Consultancy ( http://' 



'www.ffconsultancy.com I 



16 




I . 1 . 1 . 1 . 1 . 1 . 1 . 1 o 



FIG. 2: Partial structure factors Sa/siQ) of (GeSe3)o.9oAgo.io (solid lines) and (GeSe3)o.85Ago.i5 
(dotted lines) glasses. 




FIG. 3: (Color Online) Three-dimensional r — k diagrams of the instantaneous amplitude A of the 
continuous wavelet transform of the experimental structure factors^^ for (GeSe3)o.9oAgo.io (top) 
and (GeSe3)o.85Ago.i5 (bottom) glasses, using the methods of Harrop and coworkers'^ . 
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FIG. 4: Partial pair correlation functions ga/3(r) of (GeSe3)o.9oAgo.io (solid lines) and 
(GeSe3)o.85Ago.i5 (dotted lines) glasses. 
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FIG. 5: Electronic density of states and species projected electronic density of states for Se, Ge, 
and Ag for (GeSe3)o.9oAgo.io (top panel) and (GeSe3)o.85Ago.i5 (bottom panel) glasses. 
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FIG. 6: Inverse participation ratio of (GeSe3)o.9oAgo.io (top panel) and (GeSe3)o.85Ago.i5 (bottom 
panel) glasses. The vertical dot-dashed line indicates the position of the Fermi level. 
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FIG. 7: Mean square displacement for all of the atomic species in (GeSe3)o.9oAgo.io (top panel) 
and (GeSe3)o.85Ago.i5 (bottom panel) glasses simulated at T=1000 K. 



19 





FIG. 9: Histogram of RMS displacements of Ag ions in both models. Left panel, 10% Ag, right 
panel 15% Ag. 39ps of constant temperature MD at lOOOK was used to accumulate these statistics. 
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FIG. 10: (Color Online) Trajectories of the most and least mobile Ag atoms in (GeSe3)o.9oAgo 
(top panel) and (GeSe3)o.85Ago.i5 (bottom panel) glasses (T=640 K). 




FIG. 11: Local density of the most and least mobile Ag atoms as a function of time 
(GeSe3)o.9oAgo.io (top panel) and (GeSe3)o.85Ago.i5 (bottom panel) glasses. 
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FIG. 12: Average displacement of the mobile Ag ions as a function of standard deviation of the 
local density in (GeSe3)o.9oAgo.io and (GeSe3)o.85Ago.i5 glasses. The figure suggests that the more 
mobile Ag ions sample a wider range of local densities (see text). 
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